Replica treatment of non-Hermitian disordered Hamiltonians 
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We employ the fermionic and bosonic replicated nonlinear a models to treat Ginibre unitary, 
symplectic, and orthogonal ensembles of non-Hermitian random matrix Hamiltonians. Using saddle 
point approach combined with Borel resummation procedure we derive the exact large-iV results 
for microscopic density of states in all three ensembles. We also obtain tails of the density of states 
as well the two-point function for the unitary ensemble. 
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I. INTRODUCTION 



In recent years the non-Hermitian random matrix 
Hamiltonians found a broad spectrum of applications in 
physical problemsn. Their applicability ranges from dy- 
namics_pf vortex linescl and gauge theory at high baryon 
densitya to description of biological populationsQ. On the 
mathematical side, the studv_of non-Hermitian ensem- 
bles was pioneered by GinibreEl, who introduced the non- 
Hermitian counterparts of the Dyson ensembles and com- 
puted correlation functions for the unitary ensemble. A 
few decades latep, Mehta extended his results to the sym- 
plectic ensembleEl. Feinberg and ZeeQ introduced the Her- 
mitization trick and derived a number of results going be- 
yond Ginibre Gaussian ensembles. The field-theoretical 
treatmentO in the supersymmetric formulationo was. in- 
troduced by Fyodorov, Sommexs_and Khoruzhenkdl3 for 
the unitary ensemble. EfetovtiH has subsequently de- 
veloped the supersymmetric method for all three non- 
Hermitian ensembles. In this paper we take an alter- 
native route of pthe field-theoretical treatment, namely 
the replica trickE-21. This method was recently proved to 
be capable of describing the-ecpselation functions of the 
Hermitian Dyson ensemblesO-T 113. So far the validity of 
the replica approach has been restricted to the asymp- 
totic regions of the correlation functions only, while the 
supersymmetry is exact everywhere. It has, however, an 
important advantage of being applicable to much broader 
class of systems, e.g. interacting disordered systemsEj or 
Calogero-Sutherland modelscd. 

Here we show that the replica strategy works for the 
Ginibre ensembles as well. Although the treatment is 
philosophically similar to that for the Dyson ensembles, it 
requires considerable technical modifications associated 
with the intrinsic chirality of the non-Hermitian sym- 
metry classes. In particular, the Riemannian symmet- 
ric manifolds which are target spaces of the nonlinear a 
model differs from the corresponding Hermitian counter- 
parts. Rather, they are subspaces of their chiral Her- 
mitian counterparts, as may be anticipated upon Her- 
mitization. This necessitates introduction of convenient 
parameterizations, calculation of the measure and finally 
an appropriate analytical continuation. We also extend 



the method by performing all order expansion around the 
saddle points. In the n — *■ limit the resulting series may 
be summed up employing the Borel resummation proce- 
dure. As a result, we are able to obtain formally exact 
results for the density of states for all three Ginibre en- 
sembles. 

The paper is organized as follows: in Section |l| we 
collect analytic results known for Ginibre ensembles, and 
preliminaries about the replica method. In Section [n| we 
derive unitary a model using fermionic replica and obtain 
corresponding density of states, including the tails near 
the boundary of the spectrum support. Section ^ ex- 
tends the unitary a model to calculation of the two-point 
function. In Section [v| we treat the density of states for 
the symplectic ensemble using fermionic replica. In Sec- 
tion VI we treat the density of states for the orthogonal 
ensembles using bosonic replica. Finally peculiarities of 
the method and possible extensions are discussed in Sec- 
tion fvnl. 



II. PRELIMINARIES 

A. Non-Hermitian Gaussian ensembles 

Throughout, this paper we adopt Mehta's conventions 
and notationsQ. Ginibre ensemble with (3 = 1 (orthogo- 
nal), (3 = 2 (unitary), or (3 = 4 (symplectic) is defined 
as ensembles of N x N matrices H over one of the fields: 
R, C, or H. Each entry of such matrices contains (3 
real components which are i.i.d. variables drawn from 
the Gaussian distributions 



dfi(H) = (aw 



N 

nn 

8=1 i,5=l 



dH, 



(«/) 



(1) 



(a = 2,1,1 for (3 = 1,2,4). 

For a real N x N matrix H ((3 = 1), the secular 
equation, det (z — H) = 0, has real coefficients, so com- 
plex eigenvalues always appear in conjugate pairs. For a 
quaternion-real N x N matrix H ((3 = 4), one can find 
a complex 2N x 2N matrix representation (also denoted 
as H for notational simplicity) satisfying 
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H = JHJ-\ J 



Iff 
-In 



(2) 



It implies that det (z — H) = det (z — H) , and thus com- 
plex eigenvalues again appear in conjugate pairs. 
Spectral correlation functions are defined as 

R p (z u ...,z p ) = (tT5 2 ( Zl -H)---tTS 2 {z p -H)), (3) 

where (• • •) denotes averaging with dfj,(H). The macro- 
scopic density of states (DoS) of H is uniform within a 
circle of radius {fl = 1,2) or V2N (/3 = 4). Below 
we collect known results for the microscopic correlation 
functions-computed by the (skew-)orthogonal polynomial 
methodtffl and by the supersymmetry methodEJ. 

Unitary EnsemblcEIH 

I |2 00 

vR^z) = [ d\e~ x {\+\z\ 2 ) N - 1 (exact), (4a) 



T(N) 



nR\{z) = 1 

7Ti?l(z) = 1 



\z\<VN), =0 (|z|>V7V), 

„-2u 2 



(4b) 



2V2vrii 



(\z\ =VN-u, l«M«ViV), 

(4c) 



vr 2 i? 2 (z, z') = l- e-l z - z 'l 2 . (|*|, \z'\ « y/N) 
Symplectic Ensembloi0 (z = x + iy, \z\ -C y/N) 



-kR 1 {z) = 2 V 2 



dX 



-2y 2 \ 



(4d) 



(5) 



Orthogonal EnsembleQ (z = x + iy, \z\ <C V^V) 

- 2 « 2A + ^%). 



TxR l {z) = 2y I 



dX 



v/TTA 



(6) 



B. Replica method 



In the Hermitian case, one uses the identity 

5(x 



— Im = — Im— — logfa; — ie) 

7T X — It 7T CtX 

lim — Im— (x — ic) n , 
n—tO irn dx 



which in the matrix context takes the formEl 

tr6(x — H) = lim — Im— det" (a; - ie - H). (7) 

Tin dx 

The two distinct versions n > and n < are 
called fermionic and bosonic replica, correspondingly. 
The fermionic replicated generating functions (det" (a; — 
H) averaged— .with dfj,(H)) of uoitaryta, orthogonal, 
symplecticESEj, and chiral unitarjO ensembles lead to 
compact nonlinear a models on U(2n)/U(n) x U(n), 



Sp(2n)/S P (n) x Sp(n), 0(2n)/0(n) x 0(n), and U(2n), 
respectively, which are FF blocks of the Riemannian sym- 
metric suflerspaces of type AIII|AIII, BDI|CII, CII|BDI, 
and A|Ac3. After appropriate parameterization of these 
symmetric spaces one obtains n-fold compact integrals, 
which may be evaluated in the saddle point approxima- 
tion. By collecting contributions from dominant as well 
as subdominant saddle points, one obtains the asymp- 
totic expressions for the correlation functions. Bosonic 
replica has been successfully .-applied to the chiral uni- 
tary ensembleO, see also Ref.Ej. 

In the non-Hermitian case, one may employ the iden- 
tity 

8 2 {z) — —d z dz\og(zz) = lim — d z dz(zz) n , 
tt n—>0 Trn 

which in our context takes the form 
1 



tr 5 2 (z - H) 



lim — <9 z <9jdet' 



{z - H)det n (z - H^). 



(8) 



This form resembles the generating function for 
a block off-diagonal (chiral) matrix obtained upon 
Hermitizationc-l The difference is that the spectral pa- 
rameter enters off-diagonally rather than diagonally as 
in the chiral case. This distinction however does not af- 
fect the resulting saddle point manifolds, which is deter- 
mined at infinitesimal spectral parameters- Therefore, 
according to Zirnbauer's celebrated resulted, the saddle 
point manifolds are expected to be (subspaces of) U(2n), 
U(2n)/0(2n), U(2n)/Sp(n) (for fermionic replica) or 
GL(2n,C)/U(2n), U*(2n)/Sp(n), GL(2n, R)/0(2n) (for 
bosonic replica) for = 2, 4, 1, respectively, which are FF 
and BB blocks of the Riemannian symmetric superspaces 
of type A|A, AII|AI, and AI|AII. 

Our aim is to obtain the DoS and the two-point func- 
tion (cf. Eqs. (Q), (||), (||)) including exponential terms, 
using the replica strategy. We employ either fermionic 
or bosonic replica, whichever is suitable for the actual 
computation. 



III. GINIBRE UNITARY ENSEMBLE: 
DENSITY OF STATES 

Consider the fermionic replicated generating function 
(positive moment of the characteristic polynomial) of 
Ginibrc unitary ensemble defined as 

Z n (z,z)= [ dHe- tlHHf det n {z- H)det n {z- H r ), 

J C NXN 



(9) 



where dH = tt~ n Ylfj d 2 H lJ . We denote 'color' indices 
by i,j = 1, • • ■ , N and 'flavor' indices by a, b = 1, . . . , n. 
In the standard way one introduces mutually indepen- 
dent Grassmann variables ip l a , ip l a , xLj Xa> an d rewrites the 
generating function in terms of an auxiliary matrix Q a t 
(a well-known reexpression of the chiral unitary ensemble 
with flavor and color indices interchangednji^J) : 
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Z n (z,z) = 



dH [ d^dxdx exp ( //'•//'• - faizS* - H^W a - j&(*5* - H^) X l a ) 

C JVX« J 

d^dipd X dx exp(^ixlxl ~ ~ zxlxl) 

trQQ+ 



dQe 

[ dQc- tlQQt det 

JC«" 



z8 a b —Qab 

nt zS al 



N 



' z -Q~ 






-L 



'ab "° Uab 



xi 



d^d^dxdx exp ( -[ip % a xl] 

dQe- trQQl det N (zz + QQi). 



(10) 



Here d^ — J\ a i dip' l a and so forth, and dQ = 

TT"" 2 Y\ab d2 Qab- A complex matrix Q e C nxn can be 
uniquely written as (singular value decomposition) 

Q = UAV, U € U(n)/U(l) n , V € U(n), 

A = diag(Ay 2 ), A a > 0, (11) 

and a Euclidean measure dQ on C nxn is related to nor- 
malized Haar measures dU on U(n)/U(l) n and dV on 
U(n) by 

n n 

rfQ = df/dt/A(A) 2 []dA a , A(A) = J|(A a - A 6 ). (12) 

a=l a>b 

As the integrand does not depend on U and V, Z n (z, z) 
is written as an n-fold integral (even for finite N), 

oo 

Z n {z,z) = / II (dKe- x "(\ a + |z| 2 H A(A) 2 , (13) 

a=1 

up to an irrelevant constant factor that approaches unity 
in the replica limit. Notice the striking resemblance be- 
tween the DoS (Q) and the replica generating function 
( |l3| ) even at finite N. This raises the question of whether 
the replica method may be applicable beyond the asymp- 



totic analysis. We come back to this issue in Sect. VII 



A. DoS in the bulk 

Now we take the large- N limit. For \z\ < y/N, the 
saddle point equation 



1 - 



N 



= 



A+|z| 2 

has a solution A — N — \z\ 2 . Therefore 



(14) 



Z n (z,z) = e-^ N -^\ (15) 

up to .an irrelevant constant given by the Selberg 
integrals, which goes to unity as n — > 0, and 0(n 2 ) terms. 
On the other hand, for \z\ > VN, the integral is domi- 
nated by the end point A = 0: 



Z n (z,z) = \z\ 



2nN 



(16) 



Using relation (||), one obtains for the DoS 

ivR^z) = \im -d z d- z Z n {z,z) = I 1 | z j < ^ . (17) 
n—>0 n i \z\ > v N 

Therefore N complex eigenvalues are uniformly dis- 
tributed within a circle of radius ViV, in agreement with 
the exact result (4b). 



B. DoS at the edge 

If \z\ is close to the edge of the circle, such that 
\z\ = \/N — u, where 1 "^u <C VN, one may work out cor- 
rections to the uniform DoS. In this case both the saddle 
point at A a = N — \z\ 2 ~ and the end point A a = 

contribute to the partition function. Summing up all the 
contributions one finds 



Z n (z,z) 



_-n(N-\z\ 2 )' 



p=0 



( n )(-) p (2ViV U ) 2p(n - p) e- 2 ^ 2 



P ^ 11 — p 

n(dA a e-^ A »)A(A) 2 x 

a=l " 6=1 



n(^A 



e 2N 



A(A) 2 . 
(18) 



Employing Selberg integrals and neglecting constant fac- 
tors which go to unity in the n — ► limit, one obtains 



Z n (z,z) = 

-n(N-\z\ 2 ) 



(19) 



where 



I» 2 



r(n-o + l) 

a— 1 v ' 



(20) 



Since D p = for p > n, one may extend summation 
over p to infinity in Eq. (|lj). We then perform the an- 
alytic continuation n — > 0. To this end we notice that 
D p — 0(n p ) in the limit n — > 0. Therefore only the 
terms with p — and p — 1 contribute to the partition 
function in the small-n limit. As a result, one finds 



3 



Z n (z,z) = c 



_ n(iV -| z | 2 ) 



1 - n 



2tt(2m) 3 



(21) 



Employing finally Eq. (g) one obtains for the DoS close 
to the edge of the spectrum support 



ttRi(z) = lim — d z dzZ n (z, z) 
n— >0 n 



-2u* 



-2u* 



2 V2iru 



O 



(22) 



This asymptotics agree with the exact result ([k]) within 
the validity of our approximation 1 -C m 2 (<C N). The 
tail of the DoS outside the circle, \z\ > \/N can not be 
obtained in the fermionic replica. Indeed, in this case 
the saddle point is situated at negative A and the con- 
tour of integration cannot be deformed to pass through 
it. As a result the integrals are dominated by the end 
point A — 0, which leads to the zero DoS for \z\ > yf~N . 
It is possible that the bosonic replica a-, model may be 
capable to produce these tails (c.f. Ref.EZI). 



IV. GINIBRE UNITARY ENSEMBLE: 
TWO-LEVEL CORRELATION 

Consider the fermionic replicated generating function 
with two spectral parameters, 



Z n { z l, z 2, Z2) 



dHe -trHHl x 



QJV X N 



(23) 



det ™ {z x - H) det n {z x - i? + ) det " (z 2 - H) det " (z 2 - H f ) . 

Introducing the enlarged flavor indices A, B = 1, . . . , 2n 
and a 2n x 2n matrix Z, 



Z = 



ziK 

Z 2 &r 



one rewrites the generating function in terms of the aux- 
iliary matrix Qab as in the previous section, 



Z n (z u zi;z2,z 2 )= [ dQe~ U QQi ' det 



N 



Z -Q 
Q f Z 



(24) 



with dQ = 7T- 4 " 2 \[ 2 l B d 2 Q AB . 

Hereafter we concentrate on the center of the circle 
\zi |, | Z2 1 < V~N. At z\ = 2 2 = 0, the large-TV saddle 
point equation 



Q*-NQ\QQ 



o, 



is solved by 



Q = VNU, U € U(2n). 



(25) 



(26) 



For zi, Z2 finite and of order unity, one can write for such 

Q, 



det 



Z -Q 
Qt Z 



exp 



<'M> ( ( 2n\z\ 2 




where 



zi + z 2 



U) = Z\ — Z2, S — 











Therefore, 

Z n (zi, zi; z 2 , z 2 ) 



-2n(N-\z\ 2 ) 



U(2n) 



dll exp L -^- tr UsU L s 



(28) 



up to an irrelevant constant factor that approaches unity 
in the replica limit. Here dU is a Haar measure on U(2n). 
The above integrand is invariant under 



U -> U 



u 
u' 



, u,u' G U(n), 



so the saddle point manifold shrinks down to 
U(2n)/U(n) x U(n). Adopting the parameterization 
of this coset manifold employed by Verbaarschot and 
ZirnbauerEi 



UsU~ 



u 
v 



cos sin e l * 
sin e~ J * — cos 



u- 1 
V' 1 



= diag(0 a ), * = diag(^), u, v € U(n)/U(l) f 

n n 

dU = dudv T\d<j) a TT dcosOa A(cos#) 2 , 



a=l 

and substituting 



a=l 



tr[/W _1 s = 2tr cos© = 2^cos6< a , 



(29) 
(30) 

(31) 



one obtains 



Z n (zi,zi;z 2 ,z 2 ) 
l 

-2n(JV- 



z|2) /n(dA Q e^)A(A) 2 , 



(32) 



where A Q = cosO a . 



The replica limit of the above n-fold integEal-Jas 



been discussed, extensively by many authors tf3E3£3 'E2l 
(see also Ref.cJ for its asymptotic analysis). The inte- 
gral J2lj| ) over Grassmann manifold U(2n)/U(n) x U(n) 
with a height function trUsU~ 1 s satisfies the criteria, 
of the Duistermaat-Heckman localization theoren£3'EJ 
and thus the saddle point method applied to Eq. ( psj ) 
is exact for an arbitrary value of \w\ 2 . For positive large 



4 



t = \uj\ 2 /2, the integral ( |32| ) may be 'approximated' by 
contributions from the end points A a = ±1, 



The two-point function is given by 



1 



= E(-) P ( F n(l)) 
p=0 



A(A) 2 

e (n-2p)t 



7r i? 2 (zi,z 2 ) = lim —fd Zl dz 1 d Z2 dg 2 Z n (zi,zi;z 2 ,z 2 ) 



l _ e -\z 1 -z 2 \ 2 +Q 



(38) 



where 



n 



(2t)( n_ P) 2 +P 2 



r(l + a/fe) 



P/ f=ir(l + (n-a + l)/fe) 



(33) 



(34) 



for p > and F®(k) = 1. One may extend summation 
over p to infinity, since the binomial coefficient (™j = 
for p > n. After this one may perform analytical contin- 
uation, employing the small-n limit, 

P a=l r ^ ~ ^ ~ 1 )/ fe ) 

Accordingly the integral reads, for small n, 



n(dA 8 e*-) A(A) 5 



3 (ra-2)t 



a=l 



(2ty 



(2t) 



n 2 -2n+2 " 



(36) 



This agrees with the exact result ( 4d ) within the validity 
of approximation 1 <gC \zi — z 2 \ 2 {<^. N). The discon- 
nected part 1 = tt 2 Ri(zi)Ri(z 2 ) is given by the replica 

symmetric saddle point p — 0, whereas the connected 
ii 2 .. 

part — e~' zl ~ Z2 ' is given by the replica nonsymmetric 
saddle point, p = 1. In the latter case, the (n — l)-fold 
integral around A = 1 and a one-fold integral around 
A = — 1 'interact' through the Van der Monde determi- 
nant. Taking this interaction into account in the per- 
turbative manner, one finds that all terms proportional 
to e~' zl_Z2 ' l\z\ — z 2 \ 2k (k > 1) are canceled. We have 
checked this statement explicitly to the order k = 1. In 
general it is a manifestation of the Duistermaat-Heckman 
localization theorem^. 



V. GINIBRE SYMPLECTIC ENSEMBLE 

For a 2N x 2N complex matrix representation of an 
N x N quaternion-real matrix H, the constraint (|^) may 
be solved by 



Using this formula, one obtains 



Z n (zi,z 1 ;z 2 ,z 2 ) 



H = 



S_ T 
—T S 



, S,TgC 



NxN 



(39) 



-2n(N-\z\ 2 ) 



e «M 2 /2 

I, ,|2n 2 



e ("- 2 )M 2 /2 

| w |2(n2-2n+2) 



(37) 



The fermionic replicated generating function of Ginibrc 
symplectic ensemble thus reads 



Z n (z, z) 



[ dHc~? tlHHt det"(z- iJ)det n (z- H) 



dS dT e" tr ss1 - trTTt det'< 

QJV X N 

-2N 2 rr^ ^2 



S_ T 
-T S 



clef 



S_ T 
-T S 



where dH = dSdT = Tt~ 2N FJij d 2 S l H 2 T l K Once again we introduce enlarged flavor indices A, B = 1, 
express the generating function in terms of the auxiliary matrix Qab , 



Z n (z, z)= dSdT dtpd^dxdx exp -SFS* 3 - T %1 T l] - Xa ] 

7 C NXJV 



5 ij 
5 ij 



Zai 



Sal 



(40) 
, 2n and 

Xb 



dipdipdxdxexp {-i\} % A Z ab i\j b - Xa z abXb ~ ^a^aXbXb + Xa^aXb^b 
dQ / d'ipdipdxdx 



M J 

CX P \-\QabQaB + \Qab{$aXb + $bXa) + ^QbaUa^B + Xb^A.) ~ ^A Z AB^ B ~ X A Z ABXb 



dQ e 



4 trQQt 



IM 



di/jdipdxdxexp -[ip l A x\] 



>ab 



Zab 



Zab Q 



ab 



Xb 



dQe 



det 



JV 



M 



Q -Z 
Z Qt 



(41) 



5 



where z = x + iy, Z = x\.m + iys. Here dQ = 
(2ir)~ 2n ~ n Y[a>b d 2 QAB, and the integration domain 
M is a set oT 2n x 2n complex symmetric matri- 
ces, because its antisymmetric part decouples from the 
fcrmionic bilinear. A similar procedure was adopted for 
the chiral symplectic ensembleL 2 ]. 

Hereafter we concentrate on the center of the circle 
\z\ <C V2N. The large- N saddle point equation 



1 



Qt-JVQt(QQt)-l = 0) 



(42) 



is solved by 

Q = V2NU, U G U(2n). (43) 
For z finite and of order unity, one can write for such Q, 

( 1 



det 



Q -Z 
Z Qt 

1 



exp 



N 



exp 
,,2 



4A 



■tr 





uz 



-\j- l z 





y tr UsU~ x s 



(44) 



Therefore, 

Z n (z, z) = c 



n(2N-x 2 ) 



dU exp 



•tr UsU^'s 



(45) 



up to an irrelevant constant factor that approaches unity 
in the replica limit. Here dU is a Haar measure on a 
group quotient U(2n)/0(2n) that is isomorphic to the 
space of 2n x 2n unitary symmetric matrices through a 
canonical projection g i— > [/ = g.g T . Since the integrand 
is invariant under 



U 



u 
u' 



U 



u 1 
u' T 



, u,«'eU(n), (46) 



the saddle point manifold shrinks down to the intersec- 
tion of U(2n)/0(2n) and.U(2n)/U(n) x U(n). Following 
Zirnbauer and HaldaneLJ, one can parameterize the ma- 
trix U and the measure dU as 



U 



u 
u 



cos i sin 
i sin © cos © 



u 1 
« T 



© = diag(0 a ), u,v € U(ra), 
= dwdu JJ dsin9 a |A(cos 2 6»)|. 



(47) 



(48) 



Computation of the above Jacobian is detailed in Ap- 
pendix A. By substituting 



tr UsU^ 1 s = 4^ cos 2 9 a - 2n, 



(49) 



o=l 



we obtain the generating function in the form of the 71- 
fold integral 



Z n (z, z) = e-^-1,1 2 ) j f[ e -V^ |A(A)| , 



(50) 



where A a = cos 9 a . Again we have encountered a strik- 
ing resemblance between the DoS (||) and the generating 
function @. 

For y 2 ^> 1 this integral is dominated by vicinities 
of the two end points. By considering t = 2y 2 to have a 
negative imaginary part (possibly infinitesimal) , one may 
deform the integration contour to the upper half plane as 
[0,1] — > [0, ioo) + (l+ioo, 1]. Employing Selberg integrals, 
one may evaluate all the contributions with p integrals 
along [1,1 + ioo) taken in the vicinity of A = 1, whereas 
the remaining n — p integrals along [0, zoo) coming from 
the vicinity of A = 0. Whenever the end point A a = is 
adopted for some a, one may safely expand (1 — A a ) -1 / 2 
in A a , while it ought to be treated exactly when the other 
end point A a = 1 is chosen. Meanwhile we approximate 
(1 — A a ) -1 / 2 ~ 1 for the former case and find 



r " 

i n (t) = j n 



|A(A)| 



r(l±2=2) (-l)Te-f* 



E(^ P (2)) 2 



p=0 



(„-p)(n-p + l) p 2 
t 2 +2 



(51) 



As in previous sections, one may extend summation over 
p to infinity and perform analytical continuation, em- 
ploying the small-n limit (p5|). Only p = and p = 1 
terms contribute in the replica limit n — > 0, since p > 2 
terms are of the order 0(n 2 ). As a result, 



/ n (t)~t-»/2 + n !V^-l/2 e -t_ 



(52) 



In order to obtain perturbative (i.e. powers of l/t) cor- 
rections, one expand the (1 — A a ) -1 / 2 factor in powers of 
A a for the replica-symmetric (p = 0) term, 



1/2 



/ n(dA a e- tA »)|A(A)in(l-A b )- 1 / 

g a=l 6=1 

ZOO / \ 

J U(dX a e-^)\A(\)\e^(j2^y (53) 



where u>k = Y^2=i^a are elemantary symmetric polyno- 
mials. In Appendix B we prove by loop equationsEJ that 
the RHS of Eq. @ is equal to 



G 



(54) 



-1/2 



ioo / \ 

r. n / 00 k\ 

/n(rfA a e-^)|A(A)|exp £ |1 

Q 0=1 \fe=l / 

4 ~ n 

= / n( rf A e-^)|A(A)|(l-^. 

q a=l b=l 

ioo 
P n 

= (l + d t )-^j JJ(dA e- tA -)|A(A)| 



plus terms of order 0(n 2 ). Here (1 + dt) 1 I 2 stands for 
a formal series 



k=0 



(2k- 1)!! d k 
(-2) k k\ dt k ' 



(55) 



Accordingly the generating function reads 
Z n (z,z) = c-^ N ~^I n (2y 2 ), 



-l/2 e -t 



(56) 



The DoS is then given by 
irRi(z) = lim — d z d g Z n (z, z) 

n— >0 rt 

= 1- ^{dutjdtil + dty^logt + iV^e- 1 , (57) 

i = 2y 2 , up to terms of order 0(i _1 / 2 e - *). On the other 
hand, one can evaluate the asymptotic expansion of the 
exact result (|j) for Rey 2 S> 1 and Imy 2 < by using 
the same saddle point method (it is formally equivalent 
to the calculation given above for n = 1), 



irR^z) = t 



:e-* A -e-' 



v/^A 7 



-tx' 



(58) 



One may confirm the equality of Eqs. (|51]) and (j58|) either 
by using a canonical commutation relation [dt , t] = 1 to 



bring e.g. all t's to the right of cVs, or by explicit sub- 
stitution of Eq. ( [55|) to derive an identical asymptotic 
series 



■kRi(z) 



^ (2k- 



1)!! 



t W 



+ iy 2ny e 



(59) 



The formal Borel resummation of the series, Eq. (|5S|), 
leads the exact result for the DoS, Eq. (|J). Therefore 
our fermionic replica method gives correctly all orders of 
perturbative terms as well as the nonperturbative term 
of the DoS. Although the imaginary terms ±i\/2Try e~ 2y 
are artifacts of locating y 2 in the positive side of the lower 
or upper half plane and are obviously absent for real pos- 
itive y 2 (i.e. the positive real axis is the Stokes line for 
the asymptotic expansions), it is their presence and co- 
incidence that enables us to identify the two expressions: 
( pp| ) with a Laplacian applied and Q , viewed as complex 
functions of y 2 . 



VI. GINIBRE ORTHOGONAL ENSEMBLE 

For a reason that will be mentioned below, we consider 
the bosonic replicated generating function (negative mo- 
ment of the characteristic polynomial) for the case of 
Ginibre orthogonal ensemble: 

Z- n (z, z) = 

/ dHc-i tTHHT det~ n (z- H)det~ n (z- H), (60) 

dH = (2tt)- n II ' A gain we introduce enlarged 

flavor indices A, B = 1, . . . , 2n and a 2n x 2n matrix Z. 
To calculate a negative power of the determinant ( |60| ) 
with the help of commuting variables (ft, </>, one has to 
ensure convergence bii-choosing Im z > and inserting a 
matrix y/s defined asEd 



t n 

il n 



as a result: 



Z- n (z,z) 



dH / exp --H»H» + i^ A 4l AB (Z BC 8^ - S BC H^)^ CD ^> D 

R NxN J C 2nXN V Z 



d(j)d(j) exp --^aSab^b^cScd^d + i^aVs 'ab z bcVscd4>d 

Q2nx» \ Z 



dQdQe-^ tr ' 



M 



C 2„XN 



Qbc 


Zbc 




'4>h 




_ Zbc 


Qbc _ 


VsCD 


¥ _ 


) 



dQdQe~i tlQ Qdet- 



M 



Q Z 
Z Q 



(61) 
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Here 



(27T) 



~2Nn 



n 



U d2( t>A> 



dQdQ 



(2n)~ 2n ~ n riA>fl ^QabcPQab, and the integration 
domain M is a set of pairs of 2n x 2n complex symmet- 
ric matrices (Q, Q) that are mutually Hermitian conju- 
gated, Q = Q>. The restriction to symmetric matrices 
is because their antisymmetric parts decouple from the 
bosonic bilinears. 

Hereafter we concentrate on the center of the circle, 
\z\ <C \/N. The large- A'' saddle point equations 

^(QQ)- 1 Q = o (62) 



imply 



QQ 



-Nl 



2 1 1 ■ 



(63) 



Clearly this saddle point manifold lies outside M., for 
QQ would have to be positive definite. For this sad- 
dle point to be admissible, one must deform the origi- 
nal integration domain M = {Qab ^iQab 
Q ± = {Q± Q)/2, into M = {Q\ B € TL,Q AB € C AB }, 
where the contour Cab is shown in Figure 1. 




FIG. 1. Schematic view of the integration contours in the 
complex Q^b plane. The original contour (imaginary axis) 
can be deformed to Cab that passes through the saddle point 
(•) without encircling the poles (o) of the determinant, each 
of which appears to the left of the saddle point and symmet- 
rically in four quadrants. 

The width of the deformed section of the contour scales as 
1/s/N. Therefore, as long as y ^> 1/VN — > 0, the poles 



of det N / 2 g q are away from the deformed contour. 
Then the saddle point equations can be solved by 

Q = v^NU, Q = -%/iV[/-\ (64) 
C/eGL(2n,R), U = U T . 

For z small and finite of order unity, one can write for 
such Q and Q, 



det 



Q Z 
Z Q 



exp 



1 

2N 



<tr 



U^Z 
-UZ 



exp ( (2nx 2 - y 2 tr UsU~ 1 f 



(65) 



where z — x + iy, Z = x^n + iy s - 

After having deformed the contour of the integration 
over (Q,Q), one must re-ensure convergence of the <p inte- 
grals in Eq. (61 ) by keeping the real part of the quadratic 
form 



' Q Z~ 






Z Q 







y^U 2 +™Y.^«\ 2 -\¥+^) 
l 

> + Q 



A=l 

-[Re0 Im^Jv's 



Q- 





' Recj)' 




Im (j) 



(66) 



to be negative definite. A sufficient condition is 
{VsQ + Vs)ab, i(VsQ~Vs)AB e R 



~s = ±Q ± <^ Qs = sC 



(67) 



This condition is actually necessary for the real part of 
the quadratic form to be negative definite for arbitrarily 
small y > 0. At the saddle point Q = -Q^ 1 = y/NU, it 
implies a constraint UsUs = — 1. Therefore, 

Z- n (z,z) =e n ( JV - l2 > J dU exp (^ixUaU^aj , (68) 



V = {U e GL(2n, R) | U = U T , (Us) 2 



-1}, 



(69) 



up to an irrelevant constant factor that approaches unity 
in the replica limit. Here dU is a Haar measure on the 
symmetric space of symmetric GL(2n, R) matrices, con- 
sisting of 2n + 1 disconnected components with 2n — k 
positive and k negative eigenvalues. The group quotient 
GL(2n, R)/0(2n), which is the saddle point manifold of 
the bosonic replicated a model for the chiral Gaussian 
orthogonal ensemble, is isomorphic to the k = compo- 
nent by a canonical projection j nt (J = gg T , but the 
component used here as the integration domain I? is a 
subspace of the k — n component. This difference be- 
tween non-Hermitian and chiral Hermitian ensembles is 
not surprising, because the symmetry breaking parame- 
ter z plays a critical role as a regularization parameter, 
either in the bosonic replica method or in the bosonic 
part of the supersymmetry method in general. 
Eq. (|68|) is invariant under 



U 



u 
u' 



U 







u 



u,u' eO(n), (70) 



so the saddle point manifold shrinks down to the intersec- 
tion of T> and GL(2n, R)/0(n) x O(n). A generic solution 
to the constraint (p9|) is given by 



U 



u sinh© cosh© 
v cosh© sinh© 

= diag(# a ), u,v e O(n). 





" T 

u 


" 







T 
V 1 



(71) 
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Comparing Eq. ( |7l| ) with (^), the Haar measure dll in 
terms of this parametrization is formally obtained from 
Eq. by a replacement — > tt/2 — i@. Accordingly, 

n 

<2[/ = dwdv dcosh6» a |A(sinh 2 6»)|. (72) 

a=l 

By substituting 

n 

tr UsU~ 1 s = - tr £/ 2 = -4 ^ sinh 2 6 a -2n, (73) 
one finally obtains 

oo 

Z_„(z, z) = e^-W / II (^7== e" VA «) I A(A)| , 
o a=1 

(74) 

where A Q = sinh 2 a . Note that an identical expression 
also follows from the case with y = Imz < 0. Once again 
we have encountered the resemblance between the DoS 
(g) and the generating function (Q). 

For y 2 3> 1 /JV this integral is dominated by the vicinity 
of the end point A a = 0. By temporarily approximating 
(1 + A a ) -1 / 2 ~ 1, one obtains 

Z_ n (z, z) ~ e "( JV -l»l a )t- n / 2 , t = 2y 2 > 0. (75) 

In order to obtain the perturbative corrections in l/t, 
one expands the (1 + A a ) -1 / 2 factor in powers of A a as in 
Sect. [v|. Accordingly the generating function reads 

Z_ n (z, z) = e^CJV— 1*1") (i _ d t )- l 'H- n ' 2 . (76) 

The DoS is then given by 

ttRx(z) = lim d z dzZ^ n (z,z) 

n— >0 —n 

= l+ l -{d u t}dt{l~d t )- 1 ' 2 \ogt. (77) 

On the other hand, one can evaluate the asymptotic ex- 
pansion of the smooth part of the exact result for 
Rey 2 3> 1 (this calculation is formally equivalent to the 
above given one for n = 1), 

7rR 1 (z)=t(l-d t )- 1 /H- 1 (t>0). (78) 

One can confirm the equality of Eqs. ( JT7| ) and ([78]) either 
by using [<9 t , t] = 1 or by explicitly deriving an identical 
asymptotic series 

The formal Borel resummation of this series leads to the 
smooth part of the exact result (0). The S function peak 



of the DoS localized at y = cannot directly be cap- 
tured by our bosonic replica treatment that requires a 
finite imaginary part of ztQ. However it can be restored 
from the Borel-reconstructed smooth part through a nor- 
malization sum rule 

oo 

J (7rR 1 (y)-l)dy = 1 (80) 

— oo 

with a correct coefficient ^tt. Therefore the bosonic 
replica method correctly gives all orders of the pertur- 
bative terms as well as the absence of exponentially de- 
caying terms (and the Stokes line at the positive real y 2 
axis), and indirectly the 6 function term of the DoS. 

Finally we comment on the fermionic replica treatment 
of the orthogonal ensemble. Following the same steps as 
m Sect. 0, one may obtain Eq. (|45|) with the integration 
domain being 2n x 2n anti symmetric unitary matrices, 
which is isomorphic to a group quotient U(2n)/Sp(n) 
by a canonical projection g i— > U — gJg T (accompa- 
nied by a trivial change 2N — * N). At present we are 
not aware of an appropriate parameterization, such as 
Eqs. ( |47|) or (fnj), of the intersection of this coset and 
U(2n) /U(n) x U(n). Consequently these generating func- 
tions are not available in the form of n-fold integrals as 
in Eqs. @, ©, ©, and @ that are suitable for 
the asymptotic analyses. This difficulty also arises in 
the bosonic replica treatment of the symplectic ensem- 
ble, whose saddle point manifold is an intersection of 
V* (2n)/Sp(n) and U* (2n) /U* (n) X U* (n) . It is currently 
not clear to us whether this difficulty is merely a technical 
obstacle or has an essential reason behind it. 



VII. DISCUSSIONS 

We have derived the fermionic or bosonic replicated 
nonlinear a models for Ginibre unitary, symplectic, and 
orthogonal ensembles. The corresponding symmetric 
manifolds are subspaces of U(2n), U(2n)/0(2n), and of a 
space adjacent to GL(2n, R)/0(2n) correspondingly, in 
conformity with Zirnbauer's classification. The proper 
parameterization of the manifolds and subsequent inte- 
gration of irrelevant degrees of freedom lefts one with the 
n-fold compact or noncompact angular integral. Per- 
forming the asymptotic analysis of such integral, one ob- 
tains sums over variety of possible saddle points both 
replica-symmetric and replica nonsymmetric. After ana- 
lytical continuation, n — > 0, only replica-symmetric con- 
tribution and possibly few nonsymmetric one survive. 
The former are responsible for the perturbative (poly- 
nomial in inverse large parameter), while the latter for 
the nonperturbative (essentially singular) contributions 
to the correlation functions. For the DoS calculation, it is 
possible to evaluate the infinite order perturbative expan- 
sions near these saddle points. One may also show that 
the Borel resummation of such expansions after n — > 







limit leads to the exact results valid for arbitrary spectral 
parameter. 

We recapitulate our findings for the non-Hermitian en- 
sembles at hand. The bulk of the DoS in the unitary 
case (constant) is fully given by the replica-symmetric 
saddle point. The nonsymmctric contributions show up 
only if one approaches the very edge of the spectral sup- 
port, giving the exponentially small deviations from the 
abrupt behavior on the border of the circle. The irre- 
ducible part of the two-point correlation function (expo- 
nentially decaying with the distance between the eigen- 
values) is fully given by a single replica nonsymmetric 
saddle point. In the symplectic case the situation is dif- 
ferent. In addition to the symmetric contribution that 
is determined to all polynomial orders, the nonsymmet- 
ric contribution exists already on the level of the bulk 
DoS. It brings an exponentially small correction (not 
oscillatory as in the case of Dyson ensembles) to the 
symmetric contribution, which is crucial in making con- 
tact with the exact result. In the orthogonal case, there 
only exists a single replica-symmetric saddle point that 
gives all polynomial orders of the DoS, and the expo- 
nential corrections to the symmetric contribution are ab- 
sent. One can restore the S functional peak from the 
smooth part of DoS reconstructed by the Borel resumma- 
tion. The presence/ absence of exponential corrections is 
a direct consequence of compact/noncompactness of the 
saddle point manifold which is the target space of the 
fcrmionic/bosonic replicated nonlinear a models. 

Finally we list possible lines of future investigation: 

• The resemblance generally observed between replicated 
generating functions (|l3|), (|50|), ( |74| ) and corresponding 
DoS (|4a|) , (|J), (||), as well as a surprising small-n reduc- 
tion (p9| ) of the formers to essentially one-component 
integrals, makes us speculate that the replica treatment 
may be applicable to Ginibre ensembles beyond Borel 
resummation or asymptotic analyses (for large symme- 
try breaking parameters \u>\, y or at large N) and could 
indeed produce full exact results. Presumably it could 
be done by applying a complete set of loop equations to 
the n-foldj integrals to reduce them directly to one-fold 
integralsEZl. 

• Application to off-diagonal parts of the remaining two 
Altland-Zirnbauer superconducting chiral classesES, i.e. 
non-Hcrmitian complex symmetric and complex selfdual 
matricest-3, should be interesting in its own right, as an- 
alytic results have not yet been obtained via other meth- 
ods. It would require appropriate parameterization of 
the FF or BB blocks of the Riemannian symmetric su- 
perspaces D|C and C|D. 

• Starting from Hatano-Nelson disordered rjayniltonian, 



either with strong or weak non-HermiticityOO, one may 
obtain a finite-dimensional nonlinear a model of the cor- 
responding symmetry class. Our formulation and results 
should serve as a ground upon which non-Hermitian dis- 
ordered Hamiltonians in the diffusive regime (with or 
without interactions) may be treated. 
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APPENDIX A: JACOBIAN 

In this Appendix we compute a Jacobian associated 
with the parameterization (4?|) following Ref.t3, Sect. 
VIIF and App. A., where it was proven that 

(Zirnbauer-Haldane): Let G be a Lie group, G e a 
subgroup that commutes with s, K another subgroup 
whose right action on G leaves a function / of G in- 
variant, A a maximal Abelian subgroup for the Cartan 
decomposition of G with respect to K, M a subgroup of 
G e DK that commutes with all elements of A, and A + an 
open subset of A such that the map (j> : (G e /M) x A + — > 
G/K, (<?M, a) i— > gaK is bijective. Let dg be a Haar mea- 
sure on G, dgx a Haar measure on G/K, and da a Eu- 
clidean measure on A. Let T(») denote a Lie algebra for 
a Lie group and a tangent space at the origin for a group 
quotient. T(G) = T(G/K) T(K) decomposes into the 
even ( e ) and odd ( Q ) elements with respect to an involun- 
tary automorphism g t— > sgs. As the action of ad(loget), 
a e A+, maps T{K) e + T{G/K) -> T(G/K) a + T(K) e 
and T{K) + T(G/K) e -> T(G/K) e + T(K) a , the set 
of positive roots A + of G decomposes into A+ and A+ 
accordingly. Then the Jacobian J(a) associated with the 
map 4>, 

[ f(gK)dg K = f ( f f(gaK)dg) J(a)da, (Al) 

JG/K JA+ \JG a I 

is given by 

J(a) = J^J sinha(loga) j J cosha(loga). (A2) 

In the case of fermionic replica for Ginibre symplectic 
ensemble, we take 



G = U(2n), K = 0(2n), G e = U(n) x U(n), 

cos(0/2) isin(0/2) 
zsin(0/2) cos(0/2) 



( l 


' 


" 




a = CXp I 2 








)- 



= diag(0 a ), da=]Jde a , 



f(gK) = exp ( -— trgg T s(gg T ) 1 s 



(A3) 
(A4) 



and identify gg T in Eq. (A4) with U in Eq. (45). For 

9 e G e , 
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tr {ga)(ga) T s((ga)(ga) 
cos i sin 

tr 

i sin cos 



T\ — 1 2 — 2 

s = tra so s 



1 

-1 



cos© — isin® 
— isin© cos© 



1 

-1 



cos 



(A5) 



and thus the integration over G e is trivial. Without loss of generality one may restrict {9} in a cell 6\ > 62 > • ■ • > 
6*„ > 0. The eigensystem of ad(loga)A = [log a, X] = a(loga)X is given by 



Positive Root 



(0 = 1,. 



,n) 



Eigenvector 

v _ ;>(+)*(-) 



G T(G/.fs:) e + T(K) 



± 1 



A+ 



M A + (1 < a < b < n) X A 



, _ J 6% 6™ + 6<fi6£ e T(G/Jf ) e + T(if )„ 
" I 4t } 4f - e T(G/X) G + T{K) e 



where 5^ = <^a ± <^4a+n- Accordingly, 



(A6) 



j( a )da= n^n c ° s 

cos — - — cos — - — sm — 



a— 1 a— 1 



a<h 



2 2 2 2 

n 

= Yl dsin0 a Y[ (cos 2 6> a - cos 2 9 b ) = Y[ dsin6 a A (cos 2 6), 



(A7) 



a>b 



up to an irrelevant positive constant. In a generic cell, 
one would have 



J(a)da = J^J dsm6 a |A(cos 2 



(A8) 



APPENDIX B: LOOP EQUATION 

In this Appendix we derive loop equationsS for an 
(auxiliary) Laguerre orthogonal ensemble and solve them 
in the replica limit, as it is essential in-deriving asymp- 
totic series for the generating functionstll. 

Consider a set of n random non-negative numbers 
{A a } whose unnormalized joint probability distribution 
is given by 



dp(\) = ~[[d\ a e- tx "\A(\)\. 



(Bl) 



Let (■ ■ •) denote an average with respect to dp(X). An 
expectation value of a product of elemantary symmetric 
polynomials u p = J2 a A£, 

s 

(Ylwpt), p<eN, 

i=l 

is called a loop amplitude. Loop equations (a.k.a. Vi- 
rasoro constraints^) are derived from integrals of total 
derivatives 

/°°n n n / s n \ 

o=l 6=1 b V i=l c=l / 

(B2) 



(pi = 0, 1, . . .). By applying derivatives to each of the 
factors, one obtains 



Pa + l 



s ^ Po s 

(upo n ) + 9 Y.^-p ^p n ^ ) 



p=0 



■^2pi(u Po+Pi JJ u P] ) = t(u; po+1 ]\uj Pi ). (B3) 

i=i .,-(=#) i=i 



Starting from wo = n, this set of loop equations recur- 
sively determines whole loop amplitudes. First few of 
them are 



(oj 2 ) = 



n(l + n) 
2t ' 
n(l + 2n + n 2 ) 
2t 2 ' 
n(2 + in + 2n 2 + n 3 ) 



At 2 

n(8 + I9n+ 16n 2 + 5n 3 ) 



{UJ1U2) 



8t 3 

n (4 + 9n + In 2 + 3n 3 + n 4 ) 



3 _ 71(8 + 14n + 13n 2 + 9n 3 + 3n 4 + n 5 ) 
4> ~ ST 3 : 



(B4) 



Degree-counting in Eq. (B3) leads to 



(v Pl ,..., P n + 0(n 2 )) (B5) 
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for pi = 1,2, 



Now we prove a surprisingly simple formula 
v Pl ,..., Ps = \((j2n) ft = 1,2,..., (B6) 



by induction: Eq. (B6) holds at degree 1, v± = 1/2. 
Assume that it holds for all loop amplitudes of total de- 
grees J2i=oPi = An amplitude of the degree k + 1, 
u po+i,pi,—,p« w i tn + l,Pi, ■ ■ ■ ,Ps > 1 is given as a lin- 
ear combination of these amplitudes by the loop equation 
(HI. If po > 1, its LHS gives 



U PO+bPl,~ >P« 



+ 1 , Po - 1 



S 1 S 



i=0 



i=0 



thus Eq. ( ]B6| ) also holds for w po +i 
first two terms of the LHS of Eq. (IB 
in n, so one again has 



(B7) 



■j,...,p s ■ If Po — 0, the 
J) are of higher orders 



--(jXKt"-')'-^ 



(B8) 



Therefore Eq. (B6) holds for any loop amplitude of the 
total degree k + 1. End of proof. It means that the 
loop amplitude depends only on the total degree of -the 
symmetric polynomials in the small-n (replica) limited, 



(n, Pi )=(,F' p, )(i+ow) 



(B9) 



»=i 



As long as Ret > and Imt < 0, one can rotate the in- 
tegration contours of A's from [0, oo) to [0, zoo) without 
ever modifying the loop equations and violating conver- 
gence. This justifies the transition from Eq. (pa) to (54). 
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